Bistability and oscillatory motion of natural nano-membranes appearing within 

monolayer graphene on silicon dioxide 
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The recently found material graphene is a truly two-dimensional crystal and exhibits, in addi- 
tion, an extreme mechanical strength. This in combination with the high electron mobility favours 
graphene for electromechanical investigations down to the quantum limit. Here, we show that 
a monolayer of graphene on Si02 provides natural, ultra-small membranes of diameters down to 
3nm, which are caused by the intrinsic rippling of the material. Some of these nano-membranes 
can be switched hysteretically between two vertical positions using the electric field of the tip of 
a scanning tunnelling microscope (STM). They can also be forced to oscillatory motion by a low 
frequency ac-field. Using the mechanical constants determined previously [ld|, we estimate a high 
resonance frequency up to 0.4 THz. This might be favorable for quantum-electromechanics and is 
prospective for single atom mass spectrometers. 



Nano-electromechanical systems (NEMS) are promis- 
ing, e.g., as ultra-low mass detectors. One measures 
the shift in resonance frequency of a micro- or nano- 
object, when a particle adsorbs onto it. Mostly, sili- 
con nano-beams with sensitivities down to attograms are 
used P, 0, S For further improvement of sensitiv- 
ity, the oscillators have to become lighter, e.g., thinner. 
Graphene is one atomic layer thick and exhibits an ex- 
treme mechanical strength [J 0, l^li ISj which could be 
ideal. Resonators using double layer graphene have al- 
ready been demonstrated implying a mass sensitivity of 
4.5 X 10~^^ kg fl]. The highest sensitivity so far is found 
for a double wall carbon nanotube, which exhibits reso- 
nance frequency shifts compatible with the mass of one 
gold atom 

In this work, we study a graphene monolayer on Si02 
by STM. We find movable areas within the valleys of 
the intrinsic rippling fs', ^ exhibiting an extremely small 
size of 6 — 21 nm-^ (230 — 800 carbon atoms). They show 
either a hysteretic or a reversible deflection in response 
to the tip bias or a change in tip-graphene distance. The 
areas can be forced to oscillatory motion by an ac-field. 
Bistability and oscillatory motion are reproduced using 
the clamped membrane model [HI including elec- 
trostatic and van-der-Waals (vdW) forces of tip and sub- 
strate. We deduce resonance frequencies up to 400 GHz 
corresponding to a vibronic energy of 1.6 meV. Such a 
large value might provide easy cooling to the ground state 
leading to a novel access to quantum-mechanical manip- 
ulation of vibrons [isl [T^ or, in combination with the 
low membrane mass, to mass detection of single hydro- 
gen atoms. 



The morphology 
ously by STM gllj 



of grap hene has been studied previ- 
I, lid . Il7| . By comparing with the mor- 



phology of Si02, we revealed that graphene can exhibit 
an intrinsic rippling not induced by the substrate. We ar- 
gued that the rippling appears, if graphene is freely sus- 
pended between hills of the Si02 [8]. The rippling (am- 
plitude: 0.3-0.5 nm) exhibits a preferential wave length of 
15 nm at 300 K, [8j slightly increased at 4.8 K to 20 nm. 
At smaller length scales, we observe additional corru- 
gation of ±70 pm. Figure [1^ shows this corrugation at 
several tip-sample distances. At large distance (top), the 
corrugated surface and the typical carbon hexagons are 
visible. Interestingly, the apparent lattice constant dif- 
fers by 14% between oppositely curved areas (Fig. [T|d). 
This is explained by the tilting of the p^-orbitals within 
the curved surface assuming an effective j3z-length of 
0.26 ± 0.05 nm in accordance with theory [l8|. 

Decreasing the tip-graphene distance (middle, Fig. 
[1^), a small bump appears within the valley, which at 
even smaller distance (bottom) increases in height and 
diameter. With respect to the neighbouring hill, the val- 
ley is lifted by 32 pm, i.e. about half of the initial height 
difference hill-valley. Within the lifted area, the atomic 
structure changes from hexagons to bumps appearing at 
every second atom position, which has been checked by 
following up lines of atomic corrugation. The symmetry 
between A and B lattice is broken, most probably due to 
in-plane compressive stress. An explanation is sketched 
in Fig. [TJ:: the electrostatic and vdW force of the tip (Fd, 
-FvdW(tg)) lift the graphene valley until compensated by 
the restoring elastic force i^mom of the membrane and the 



vdW force F, 



vdW(sg) 



of the substrate. Since the two tip 
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forces change with lateral tip position, a dynamic image 
of the lifting results as indicated by white dots in Fig. [T};. 
While a hill appears within the STM image, the mem- 
brane still maintains its valley-like shape, but with re- 
duced curvature. The resulting compressive, lateral force 
within the lifted area can be reduced by a vertical zig- 
zag atomic arrangement straightforwardly explaining the 
observed symmetry breaking between A and B lattice. 
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FIG. 1: Reversible lifting of the graphene nano- membrane, a) 3D representations of an atomically resolved monolayer graphene 
at different tip-sample distance as indicated on the right. At a large distance {U = IV, I — 0.1 nA) carbon hexagons are 
visible. Decreasing the distance, a hill appears in the centre of the valley and the hexagons are transformed to bumps at every 
second atom position, b) High-pass filtered STM image ([/ = IV and I = 0.1 nA). Circles denote the hill and valley areas 
marked in a. The apparent lattice constants show a difference of Aa — 0.03 nm due to the tilting of the vr-orbitals. c) Model 
explaining the mechanical behaviour of the nano-membrane. Due to the locally applied force between tip and graphene, the 
valley is continuously lifted while scanning the tip and the STM image shows a hill in the centre of the valley (dotted line). 



While the membrane's lifting in Fig. [T] is reversible, 
other valleys exhibit hysteresis. FigureOa shows the I{z)- 
curve of a hysteretic area at a sample voltage of J7 = IV. 
While approaching or retracting the tip, a jump in cur- 
rent / by three orders of magnitude is observed with a 
hysteresis of 50 pm. Such hysteretic membranes can be 
identified directly within constant-current images, if / is 
chosen within the bistable range. This is demonstrated 
in Fig. where a valley (hill) showing atomic hexagons 
is visible at low (high) /, while a noisy area appears in 
between. The noisy behaviour is caused by the feed- 
back loop, i.e. the tip approaches which increases the tip 
forces and induces a snapping of the membrane towards 
the tip; this increases / by three orders of magnitude and 
the tip retracts again leading to reduced tip forces and a 
snapping back of the membrane. We observed bistability 
only in valleys and never on hills. Although analysing 
several tens of valleys, we could not find any correlation 



between hysteretic or reversible behaviour and depth, 
width or curvature of the valleys. Thus, we conclude 
that graphene-substrate interactions not visible by STM 
are crucial. 

In order to model the observed behaviour, we analyse 
the involved interaction potentials. Besides the electro- 
static potential $oi between tip and graphene, the vdW 
potentials $vdw for tip/graphene and graphene/Si02 
[l9|, and the elastic potential of the membrane $mom 
are considered (see methods). They are plotted in Fig. 
[2]; as a function of vertical graphene position for a fixed 
tip-substrate distance (1.63 nm). The summed up poten- 
tial $sum exhibits two local minima representing the ob- 
served bistability. Fig. shows two nearly degenerate 
minima at U — OY and tip-substrate distance 1.46 nm 
(top). The lower minimum transforms into a saddle at 
U = 1.3 V (middle) forcing the membrane to the up- 
per minimum, i.e. 0.34 nm closer to the tip. Increasing 




FIG. 2: Hysteretic nano-membrane. a) Constant-current images of a graphene area taken at different tip-surface distances. 
At intermediate distance (/ = 3nA, 15 nA), a part of the valley exhibits unstable behaviour visible as spikes. Higher current 
leads to stable imaging again, but the valley has turned into a hill, b) I(z)-spectrum revealing hysteretic behaviour, c) 
Calculated potential energy "I>sum as a function of substrate-graphene distance Zsg at fixed tip-substrate distance. $sum consists 
of electrostatic potential $ci, van-der-Waals potentials between tip and graphene 'l"vdw(tg) and Si02 and graphene 4>vdw(sg), 
and the elastic membrane potential <E>mcm. The two local minima in ^sum indicate two metastable positions. The definition of 
distances is sketched at the top. d) Using an initial tip-substrate distance 2tg -I- Zsg = 1.46 nm, (top) the hysteretic switching in 
constant-current mode can be reproduced by first applying a voltage moving the membrane towards the tip (middle) and then 
retracting the tip to reduce current again, thereby switching the membrane back (bottom). 



now the tip-graphene distance Ztg by 0.34 nm (reduced /) 
shakes the potential back and the membrane flips back to 
its original minimum (bottom). Note that Ztg ~ 0.56 nm 
is estimated by extrapolation of I{ztg) towards the con- 
tact conduction Gq — 2e^//i [2l| leaving only the initial 
substrate-graphene distance Zgg — 0.9 nm as a fit param- 
eter. By decreasing the initial ^sg to 0.63 nm, we still 
calculate two minima, but we are not able to switch the 
membrane to the upper position by reasonable U. This 
explains the reversible behaviour observed in Fig. [TJ 

To explore the electric response of the membrane, we 
varied U. First, we measured I{U) on a hysteretic val- 
ley in comparison with Io{U) on a hill (Fig. [3^). As- 
suming the same electronic structure at both positions, 
the ratio of the two curves is given by I{U)/Io{U) = 
exp(— 2k(C/) • z+([/)), where k{U) is the electron's de- 
cay constant determined within the supplement. Thus, 
the logarithmic plot of I{U) / I{Uq) directly displays the 
deflection z+ of the valley with respect to the hill (up- 



per inset. Fig. [3^). The minimum of z+(C/) marks the 
absence of electrostatic force at the contact potential 
Uc = —0.45 V implying a tip work function of 5.11 eV 
(graphene: 4.66 eV 20]). 

Next, we apply an ac-voltage U {t) = ?7dc+C^o cos(27rz/t) 
with varying amplitude Uq at a dc-offset Udc = —0.3 V 
~ Uc and fixed tip-substrate distance. The tip is placed 
above the nano-membrane marked by the dashed circle in 
the inset of Fig. We choose /stab = 0.2 nA low enough 
to avoid snapping of the membrane. Figure [5}: shows the 
resulting I{t) and the apphed C/,„od(0 = U{t) — Udc for 
the membrane (bottom) and on a reference (hill) position 
(top). At reference, I{t) is dominated by the capacitive 
crosstalk (phase shift of 90°). At the membrane, an addi- 
tional in-phase signal, non-linear with respect to f7mod(i), 
indicates the reversible membrane movement. Fig. 
displays the corresponding <i>sum(2:sg) at the two extrema 
of electrostatic force (ztg = 0.56 nm). The membrane 
does not switch between the two minima, but oscillates 
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FIG. 3: Oscillatory deflection of the graphene nano-membrane. a) /(l7)-curve {Ustah = —0.3 V, /stab = 100 pA) recorded on a 
valley (red) and on a hill position (black) as marked in the lower inset. Upper Inset: logarithmic plot of the ratio of the two 
I{U) curves being proportional to the deflection z+(f/). The minimum marks the contact potential Uc; blue line: parabolic fit. 
b) Calculated potential ^sum for an initial substrate-graphene distance Zsg ~ 0.63 nm and tip-graphene distance Ztg = 0.56 nm. 
Even at f/ = 1.85 V the system cannot move to the right minimum, but exhibits only a reversible shift of 30 pm. c) Current 
response (black: hill, red: valley) to an oscillating tip voltage Umod (blue) measured at reference (hill) and on nano-membrane 
(valley), d) Upper inset: in-phase lock-in signal recorded on reference (black) and nano-membrane (red). Main: ratio of the two 
lock-in signals tli from the upper inset (symbols) in comparison with a fit using the model of a prestrained clamped membrane 
(blue line) [l^]. Right and upper scale show the deduced deflection amplitude Az (right axis) and the applied electric field 
amplitude _Eo (upper axis). Lower inset: STM image with the two measurement positions (white circles) and area of the 
membrane (dashed circle) marked. 



reversibly within one valley in agreement with experi- 
ment. 

Finally, we use lock-in technique to measure the de- 
flection amplitude Az with respect to C/q- The inset of 
Fig. [5Ji shows the lock-in output measured on the nano- 
membrane (red) and on a reference position (black). The 
stronger signal at the membrane again indicates its os- 
cillatory motion. To get rid of unknown parameters we 
use the ratio of the two lock-in signals 

/;+^/(Ocos(a;t)d^ 

^LI = 7f+T ■ (1) 

Io{t) cos{ujt) dt 

From tli we deduce Az numerically (see supplement) as 
displayed on the right scale of Fig. (SJi. The maximum 
amplitude is higher than expected from the model in Fig. 
This is likely caused by the fact that also the hills 



are simultaneously lifted by the tip forces, which con- 
tinuously decreases 3>vdW(sg) with respect to our model. 
Nevertheless, the displayed Az and the corresponding 
electric field amplitude Eq deduced from Uq and Az (see 
supplement) give an impression of the amplitude-field re- 
lation achievable by external excitation. 

The basic resonance frequency of the nano-membrane 
can be roughly estimated by clamped membrane the- 
ory using E^^ = 340 N/m [10|], membrane radius r = 
2.58 nm, and mass mo = 1.59 x 10~^^kg and assuming 
HOPG values for thickness h = 0.335 nm, and Poisson- 
ratio s = 0.16 [2^ : 

More adequate molecular dynamics simulations (MDS) 
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for monolayer graphene reveal = 400 GHz for a smaller 
area (4.3 nm^) exhibiting 75 % deviation from eq. [2] _23]. 
Anyhow, the nano-membranes consisting of only 200-800 
atoms are ideal resonators for ultra sensitive mass detec- 
tion. Adsorbing one hydrogen atom of mass ma would 
lead to a relative frequency-shift of 



mo/ 



-0.5 



10" 



(3) 



MDS finds quality factors Q = 2500 (273.000) for un- 
supported graphene monolayers at 300 (3) K [23| . Exper- 
iments reveal lower values for multilayer graphene prob- 
ably due to interlayer friction 1, 4, &, .25], but Q = 4000 
has been achieved using oxidized multilayer graphene at 
300 K [13 . This might be sufficient to detect frequency 
shifts of 10"'*. Measuring membrane oscillations would 
be eased by the nonlinearity of J(z) requiring only dc 
detection, but the excitation in the range of 100 GHz re- 
mains a technical challenge. 



I. METHODS 



two z^-curves with origins separated by 0.12 nm in or- 
der to describe the fact that the membrane is first later- 
ally compressed by lifting with respect to the hills, but 
is relaxed again after being higher than the surrounding 
hills. We do not include a pretension term pi3], since 
pretension is given explicitly by the vdW forces. The 
electrostatic energy $ci = 1/2 ■ C^nU^s is calculated us- 
ing Cm as the tip-nano-membrane capacitance modelled 
as a sphere in front of a plane and UeS as the voltage 
drop between tip and graphene reduced with respect to 
U due to graphene's finite carrier density. 
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III. SUPPLEMENTARY INFORMATION 

A. Determination of the initial tip/graphene 
distance z^g and the decay constant k,{U) 



The preparation of the graphene sample is done by 
mechanical exfoliation on a Si02 substrate as described 
elsewhere 0, ^26^ ■ A graphene flake containing a mono- 
layer region is identified by an optical microscope. In 
addition, the film thickness is confirmed by Raman spec- 
troscopy [sj. A gold contact surrounding the graphene 
is produced by e-beam lithography. In order to remove 
the residual resist and adsorbates as water, the sam- 
ple is heated to 170°C for a few hours, first in air, and 
then, directly before the measurement, in ultra- high vac- 
uum {p — 2 X 10""'^'' mbar). The monolayer region of 
(18 X 26) /im^ is positioned below the tip of the STM by a 
piezo motor using the control by an optical long distance 
microscope with a resolution of 5 /^m. The measurements 
are performed with a high resolution STM operating at 
T = 4.8 K [23. The tip is prepared prior to measure- 
ments by applying voltage pulses and field emission on 
the gold contact region. All STM images are measured 
in constant-current mode with the voltage U applied to 
the sample. All spectroscopic curves are measured with 
the feedback switched off at stabilization voltage J/stab 
and stabilization current /stab- Afterwards, either U or 
the tip sample distance z is changed, while measuring the 
tunnelling current /. 

Calculations of the interaction potentials acting on the 
graphene membrane are described in detail in the sup- 
plement. In short, we use a dielectric plane for Si02, 
a metallic circle of cosine shape in vertical direction for 
the membrane and a metallic W tip of parabolic shape 
(central radius: 0.7 nm) to calculate the vdW potentials 
[l9t . The elastic potential i>mem is modelled by a clamped 
membrane implying 'I'mem = A-z^ with A taken from ex- 
periment pl| and z+ being the membrane deflection. We 
assume two minima (two relaxed membrane positions) of 



The absolute distance between tip and graphene ob- 
tained after stabilization, z^g, could be determined, in 
principle, by decreasing the distance between tip and 
sample until the conductance reaches the conductance 
quantum Go = 2e^/h [2l|. Unfortunately, this does not 
work for graphene, because the graphene is lifted during 
the tip approach, even on the hills (reference positions) as 
indicated from significantly too large decay constants k 
extracted from /(z)-spectra, which have been measured 
systematically on the graphene surface. The only way 
to reasonably estimate is to calculate k — 23.3 nm~^ 
with the help of equation[5] (below) and, then, to extrapo- 
late the exponential behaviour of the conductance with 
respect to the distance by using the stabilization param- 
eters Us = -0.6 V and 1^ = 200 pA: 



Us 



2e' 



exp{-2«;(C/s)z{'g} 



^ z. 



2k{Us) 



■In 



Go Us 



0.56 nm. 



(4) 



To determine the decay constant k{U), we use a planar 
tunnelling junction with a correction factor ^ (28j : 



/2mo 



eU 



(5) 



The effective work function 0cff = (0gr + 0t)/2 = 4.89 eV 
has been calculated using the known graphene work func- 
tion of 0gr = 4.66 eV 20] and a tip work function of 
(pt = 5.11 eV derived from the measured contact poten- 
tial as displayed in Fig. 3a of the main text. The correc- 
tion factor ^ = 1.06 has been determined by fitting I{z)- 
curves measured with the same microtip on Au(lll) by 
equation O (work function: 0au = 5.31 eV ]29] ) . 
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Calculation of the interaction potentials acting 
on the nano-membrane 



In order to describe the observed behaviour of the 
nano-membrane, the involved interaction potentials are 
analyzed in detail. Besides the electrostatic potential $ei 
induced by the tip, the Casimir/van-der-Waals potentials 
<i>vdw induced by the tip and the Si02-substrate as well 
as the elastic restoring force of the membrane itself ^mcm 
are considered. 



1. Casimir/van-der-Waals potential induced by the tip 

The description of the van-der-Waals and Casimir po- 
tential <i>yriw per unit area A between two materials is 
given by [30l |: 



vdW 



4^ 



/ k^dk^ / In [1 - vrwe-^^^'s] diu (6) 



where r^^ and denote the frequency dependent reflec- 
tion coefficients of graphene and tungsten, respectively, 
fc_L is the wave number parallel to the surface, Ztg is the 
variable distance between graphene and the tip apex and 
u) is the frequency {h: Planck's constant). The reflection 
coefficient of graphene can be calculated by [13] : 



c^qVt 



with O = 6.75 X 10^m~\ c = 3 • 10^ m/s and 



(7) 



(8) 



To determine the reflection coefficient of the tungsten 
tip, the frequency dependent dielectric constant e(w) has 
to be used: 



with 



s{uj)q — k 
e{uj)q -{- k 



(9) 



(10) 



The dielectric function can be approximated knowing 
the plasma frequency of tungsten lv^^ — 9.74 x 10^^ Hz 



e{uj) = 1 + 



■^p,W 



(11) 



The total interaction potential <i>vdw(-2^) is determined by 
an integration of equation [S] over the circular area of the 
nano-membrane: 

<i>vdw(2tg(r' = 0)) = 27r^' ^I^{^z,^[r')y dr' , (12) 



tungsten tip 



graphene 
membrane 




Si02 / 



SI: Definition of distances used in the calculation of the 
interaction potentials. 



where ztg^{r') denotes the vertical distance between tip 
and graphene at a lateral distance r' measured from the 
centre of the membrane and ztg(r' = 0) indicates that 
we use the tip-graphene distance in the centre of the 
membrane as the variable for <&vdw • Assuming a circular 
membrane with the measured radius of r = 2.58 nm and 
a 2D-cosine shaped corrugation (see figure SI), as well 
as a parabolic tip with central radius R, ztg(r') can be 
described as: 



^tg(r') = — + [ztg(O) 



z'^^] cos' 



TTV 

"27 



hill 



(13) 



The vertical distance between tip apex and the hills 
surrounding the membrane zjg" is determined by analyz- 
ing the measured valley depth with respect to the sur- 
rounding hills. For the tip radius R we can only give 
an upper limit of i? = 2.3 nm determined by analyzing 
STM images of atomically resolved valleys, which would 
not have been resolved by larger tips due to convolution 
effects. The smallest possible tip radius of i? = 0.3 nm is 
given by a tetraedric alignment of the first four atoms. 
For the calculation, we use a value of i? = 0.7 nm, deter- 
mined as described below. 



2. Casimir/van-der-Waals potential induced by the 
Si02-substrate 

The interaction potential between graphene and the 
amorphous Si02-subtrate <&vdw(sg) is calculated similarly 
using equation [6V-fT0l and [T2l but replacing Ztg{r') by the 
distance between graphene and the Si02 substrate Zsg(r') 
as well by the reflection coeflicient of Si02 rsi02- 

In case of an insulating material like Si02, ^(w) is given 
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by another expression [30|: 



4- Electrostatic potential induced by the tip 



£(0) ~ 1 
1 + ^ 



(14) 



where 



1.05 X 10^^ s ^ is the main electronic ab- 



sorption frequency being within the ultra violet region 
[Ulll]. For amorphous Si02, we use the known dielec- 
tric constant at zero frequency of £(0) — 3.9 31]. To 
describe Zsg(r'), assuming a plane substrate and again a 
2D-cosine shaped membrane as sketched in Fig. SI, we 
use 



(r') 



^sgV 



^1 <-) 



with z^j" being the distance between substrate and the 
borders of the membrane (surrounding hills) and 2sg(0) 
being the distance between the substrate and the lowest 
point of the membrane. The value of Zsg(O) can vary with 
the applied voltage U or with the tip-substrate distance. 
However, the change of Zsg(O) with U or tip-substrate dis- 
tance can be measured via the tunnelling current. Only 
the offset of the initial tip-substrate distance, z^g, with- 
out electric field has to be taken as a fit parameter. 



3. Elastic membrane potential 

In order to determine the elastic potential $mcm of the 
nano-membrane, we assume a cubic force dependence of 
the deflection as given by the classical clamped mem- 
brane theory according to Komaragiri et al. [ll[ [l3| ■ The 
nanomembrane gets laterally compressed if lifted until 
the centre of the membrane is at the same height as the 
surrounding hills (see Fig. Ic of the main text). If lifted 
further the membrane gets relaxed again up to a second 
stable position above the surrounding hills. In order to 
model this behaviour, we compose the membrane poten- 
tial $mcm of two parts $+(Zsg(0)) and <i>~(zsg(0)) with 
minima vertically symmetric with respect to the position 
of largest compression leading to: 



cl>+(z. 



-0.265£;^ 



(^sg 




d\4 
21 


, if Zsg 


^ ^sg 








(■^sg 


^hill 

^sg 


d\A 
2> 


, if ^^sg 


^ ^sg 









«>"(zsg) ^ -0.265£;2 



where E^^ denotes the two dimensional Young's mod- 
ulus due to compressing or stretching of the atomic 
bonds, which has been measured previously to be E^^ = 
340 N/m [l^. The distance between the two potential 
minima d is the second free parameter of our model, only 
limited to about twice the valley depth. We found that 
d = 0.12 nm is able to reproduce the behaviour of the 
membranes displayed in Fig. 1 and 2 of the main text. 
This number is lower than the valley depth appearing in 
the STM images implying that the van-der-Waals force 
of the Si02 substrate increases the corrugation strength 
within the graphene flake. 



The electrostatic potential $ei induced by the tip is 
given by 



(16) 



with the capacitance Cm and the effective gap voltage 
Uos to be determined. We approximate the tip-sample 
system by a capacitor consisting of a sphere with radius 
R above a circular plate of radius r corresponding to 
tip and membrane, respectively. The capacitance C for 
r = oo can be calculated analytically resulting in [36| : 



C = 47re£o^sinh 

oo 

•^sinh-1 



i?cosh 



ztg + R 
R 



nR cosh f *^ ^ ^ 
R 



(17) 



with the distance between the tip apex and the infinite 
plate Ztg = ztg(O). The capacitance determined by equa- 
tion [T7] has to be modified because of the finite area of 
the nano-membrane below the tip. Therefore, we assume 
a Gaussian shape of the total charge density of the infi- 
nite plate with a maximum at r' = to be determined 
below and calculate the capacitance of the finite mem- 
brane by integration only over the circular area of the 
nano-membrane. We end up with a capacitance Cm of 



Cm — C 



1 — exp 



Cztg 



(18) 



Because of the finite charge carrier concentration of 
graphene, the effective voltage J7off between the tip and 
graphene is smaller than the applied bias- voltage U. The 
remaining voltage leads to a considerable Fermi level shift 
within the graphene until the charge carrier density is 
high enough to screen the electric field. As a consequence, 
there is a potential drop between the graphene just below 
the STM tip and the gold electrode connected to the ex- 
ternal power supply. Due to the linear dispersion relation 
of graphene the two dimensional charge carrier density n 
can be written as 13411: 



(19) 



where U := U — Uc with Uc being the contact potential 
determined in Fig. 3a of the main text, vp = 1.1 x 10^ m/s 
is the Fermi-velocity of graphene and e is the electron's 
charge. In equilibrium, the electrons screen the electric 
field E and the resulting 2D charge density can be most 
easily approximated by a plate capacitor leading to: 



ne 



eeoE{Ues) = seq 



cff 



(20) 



Thereby, z^ denotes the distance between the plate and 



the centre of mass of the lower half sphere approximating 
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the tip, which has been chosen to map the model of two 
paraUel plates to the model of a sphere and a plate. With 
the help of equations[T9]and[2Ql the effective voltage drop 
between tip and graphene becomes: 



U,ff = U 



a 



4(4)2 



aU 



(21) 



electric field and Ucs{t) is the part of the corrected volt- 
age U^{t) dropping between membrane and tip. Note, 
that the expression for is assumed not to be present 
on the reference position. The quadratic dependence of 
the tunnelling current with respect to the bias voltage 
is derived from the linear dispersion relation of graphene 
and has been checked by according fits to the measured 
/(C/)-spectra. In Fig. 3d of the main text, we display the 
lock-in ratio: 



The value of the dielectric constant for graphene on 
Si02 has been calculated previously using the image po- 
tential method and amounts to e = 2.5 [35|- The result- 
ing n calculated straightforwardly by inserting C/efi into 
equation [19] has been used self-consistently as the maxi- 
mum of the Gaussian charge density required to calculate 

Cm- 



C. Excitation of the nano- membrane by ac- voltage 

We applied an ac-voltage U{t) = Udc + Uq cos(27rz^i) 
with a varying amplitude [/q, at a dc-offset Udc and a fre- 
quency of = 1.4 kHz. In addition, we define the voltage 
W^i^t) — Udc + Uocos{2TTi't) — Uc including the contact 
potential Uc determined from Fig. 3a of the main text. 
The time dependent tunnelling currents Io(t), measured 
at the stable reference position and I{t) measured above 
the nano-membrane can be described using the linear 
graphene density of states as: 



^o(i)oc^^exp{-2K(;7(i))4} 



\U{t)\ 



(22) 



and 



U{t) 



I{t) « ^exp{-2«(t/(i)) [zj'g - z+iUcf^m } , (23) 

where z+([/eff) denotes the deflection of the nano- 
membrane with respect to its position in the absence of 



//+'^ lit) cos{2TTvt) dt 
It^'^ Io{t) cos{2Tri^t) dt 



(24) 



which can be directly used to determine the deflection 
amplitude Az numerically as displayed on the right of 
Fig. 3d of the main text. Accordingly, the upper scale 
of Fig. 3d of the main text shows the electric field am- 
plitude £'0 given by Eq = C/off,max/2tg ,min with 2:tg mill 
being the minimal distance between graphene and tip 
and C/eff,maxbeing the maximal effective voltage during 
the oscillation. 



Finally, we describe our estimate of the tip radius 
R. Therefore, we use the simplified clamped membrane 
model for the force-deflection curve, consisting of a lin- 
ear term caused by so-called pretension and a cubic term 
describing the compression of the atomic bonds given by 
E^^ 10]. The equilibrium between the electrostatic force 
and the elastic membrane force is then given by: 



2 dztg 



(Ucsit))' 



2D 



TTZ+ - E'^^im 



(25) 



where denotes the two dimensional pretension. After 
solving towards z+(?7off) numerically, we fitted the mea- 
sured riA taking and the tip radius R as the only 
free parameters. The excellent fit displayed in Fig. 3d of 
the main text results in a tip radius of i? = 0.7 nm and a 
pretension of CTq^ = 0.62 N/m. 
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